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1.0  Introduction 


Temperature  is  a  primary  variable  in  combustion  pro¬ 
cesses  and  its  accurate  measurement  is  an  essential  pre¬ 
requisite  for  a  proper  analysis  of  the  combustion  process. 
Presently  available  electrical  techniques  such  as  thermo¬ 
couples  and  resistance  thermometers  are  limited  in  their 
transient  response  and  indicate  temperatures  only  at  the 
thermometer  surface.  Optical  techniques  have  excellent 
transient  response  but  cannot  measure  temperature  in  the 
interior  of  optically  thick  carbon-rich  flames  or  combus¬ 
tion  products. 

This  report  summarizes  an  effort  to  demonstrate  the 
feasibility  of  measuring  the  temperature  of  hot  gases  by 
sensing  the  thermal  noise  power  emitted  by  the  gas. 

Thermal  noise  or  Johnson  noise  is  the  noise  produced 
by  thermal  agitation  of  charges  in  a  conductor.  The  avail 
able  thermal  noise  power  produced  in  a  resistance  is  inde¬ 
pendent  of  the  nature  and  the  value  of  the  resistance  but 
is  proportional  to  the  absolute  temperature  and  the  fre¬ 
quency  bandwidth  over  which  the  noise  is  measured. 

Utilization  of  this  universal  property  of  conductors 
leads  directly  to  a  means  of  non-intrusively  obtaining 
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transient  temperature  measurements  in  the  interior  of  a 
combustion  process.  The  measurements  are  thensodynasii- 
cally  correct  even  in  the  presence  of  particulate  siatter, 
flow  velocity  and  high  pressures. 

Noise  thermometry  in  solid  conductors  is  an  estab¬ 
lished  measurement  technique.  It  is  used  in  low  tempera¬ 
ture  research  to  measure  temperatures  in  the  millikelvin 
range  to  a  high  degree  of  accuracy  (References  1  and  2)  . 

The  nuclear  industry  employs  noise  thermometry  to  measure 
liquid  metal  temperatures  to  within  a  kelvin  by  immersing 
a  conductor  in  the  molten  flow  (References  2  and  3) .  This 
technique  has  been  considered  for  application  in  fossil 
fuel  plant  process  control  (References  4  and  5) .  These 
measurements  are  usually  made  at  low  frequencies  (<100  kHz) 
and  are  limited  in  their  time  response  because  of  the 
thermal  inertia  of  the  probes. 

The  Soviets  report  direct  measurements  of  acetylene- 
oxygen  combustion  flows  made  in  the  5  to  30  MHz  range  and 
conclude  that  high  temperature  gasdynamic  measurements  are 
feasible  (Reference  6) .  They  report  a  frequency  dependence 
of  noise  temperature  measurements  which  was  not  observed  in 
this  work. 
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2.0  Theory  of  Thermal  Noise 


Thermal  noise  emitted  from  resistive  components  has 
long  been  the  subject  of  research  since  it  is  the  limiting 
factor  for  radio  receiver  sensitivity.  The  classic  work 
by  Lawson  and  Uhlenbeck  (Reference  7)  on  receiver  noise 
for  World  War  II  radars  provides  a  clear  and  concise  theo¬ 
retical  basis  for  noise  analysis  which  is  now  an  important 
engineering  discipline. 


2.1  Derivation  of  Thermal  Noise  Spectrum 

The  noise  spectral  density  is  a  universal  function 
of  the  resistance,  temperature  and  frequency  bandwidth 
and  is  derived  from  the  fundamental  laws  of  thermodynamics. 

Consider  a  conductor  of  resistance  R  at  temperature 
T.  Because  of  the  random  motion  of  electrons  there  will 
be  small  fluctuations  of  the  voltage  across  the  ends  of 
the  conductor;  the  average  value  of  the  fluctuations  being 
zero.  Nyquist  (Reference  8)  showed  that  the  spectrum  of 
these  voltage  fluctuations  is  constant  up  to  very  high 
frequencies  (at  least  up  to  10*  MHz)  so  that  for  all  prac¬ 
tical  purposes  it  is  a  white  spectrum  with  spectral  density 
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given  by  the'  formula 

=  G(f)Af 

=  4kTRAf  (1) 

where  R  =  resistance 

Vn  =  r.m.s.  voltage 
T  =  temperature 
Af  =  frequency  bandwidth 
k  =  Boltzmann's  constant 

The  proofs  of  this  result  are  all  based  on  the  gen¬ 
eral  principles  of  statistical  mechanics  and  especially 
on  the  theorem  of  equipartition  of  energy. 

Nyquist  demonstrates  the  result  for  power  transfer 
between  two  equal  resistances  at  the  same  temperature  by 
using  transmission  line  theory  and  the  equipartition  of 
energy  (Reference  7)  . 

In  a  second  approach  a  resistance  R  possessing  a 
fluctuating  electromotive  force  E(t)  is  connected  to  an 
ideal  network  made  up  of  inductances  and  capacitances. 

By  applying  the  equipartition  theorem  to  the  electric 
and  magnetic  energy  of  the  network  the  result  of  equation 
(1)  can  be  drived  and  in  addition  it  can  be  proved  that 
E(t)  must  be  a  Gaussian  random  process  (Reference  7). 


Although  the  fundamental  result  of  equation  (1)  is 
independent  of  the  mechanism  of  conduction  it  has  also 
been  derived  from  a  specific  model  for  the  electrical 
conduction  through  a  metal. 


The  equivalent  circuit  of  a  resistance  is  a  thermal 
noise  source  of  fluctuating  voltage  Vj^  in  series  with  an 
ideal  noise  free  resistor  R.  For  a  given  amplifier  input 
resistance  (assumed  for  now  to  be  noise  free)  the  ther¬ 
mal  noise  power  appearing  at  the  amplifier  input  is 


Pi 


Pi 

4kTAf 


RRi 

(R+Ri) ^ 


(2) 


Maximum  input  power  occurs  at  matched  impedance  conditions 
(R=Ri) 


Pmax  ~  kTAf 


(3) 


and  depends  only  on  temperature  and  frequency  bandwidth. 
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2.2  Generalization  of  Theory 


The  theorem  of  Nyquist  states  the  result  of  equation 
(1)  in  more  general  form  by  considering  a  two-pole  linear 
passive  network  at  temperature  T.  If  the  impedance  between 
the  poles  is  Z  =  R(f)  +  iX(f)  where  R(f)  and  X(f)  are  in 
general  frequency  dependent,  then  it  is  shown  that 

=  4kTR(f)Af  (4) 

That  is,  the  fluctuating  emf  emitted  by  linear  passive 
network  at  temperature  T  depends  only  on  the  real  part  of 
the  network  impedance . 

The  theorem  of  Williams  shows  that  equation  (4)  can 
be  derived  by  associating  with  each  of  the  constant  resis¬ 
tances  in  the  network  a  fluctuating  emf  having  the  constant 
spectrum  of  equation  (1) .  The  result  is 

=  4k4f|zn  ?  (5) 

This  can  be  used  to  predict  what  will  happen  when  the  dif¬ 
ferent  resistances  are  no  longer  at  the  same  temperature. 


2.3  Experimental  Foundations  of  the  Theory 

Many  experimental  investigations  have  confirmed  the 
basic  Nyquist  formula  of  equation  (1).  Wilbur  (Reference 
9)  amplified  voltage  fluctuations  across  wire  wound  re¬ 
sistors  of  resistances  up  to  about  2  megohms  and  compared 
them  with  fluctuations  due  to  pure  "shot"  noise  which  is 
known  to  obey  precisely  the  theoretical  formula 


i^  =  2eIAf  (6) 

where  e  =  electronic  charge 
I  =  current 

This  work  verified  the  linear  dependence  of  G(f)  on 
R  and  T  from  liquid  air  temperature  to  380°K  and  can  al¬ 
most  be  considered  as  a  precision  measurement  of  Boltzmann's 
constant  k. 

The  most  detailed  verifications  of  equations  (4)  and 

(5)  were  made  by  Williams  (Reference  10).  In  one  series 

of  experiments  the  voltage  fluctuations  across  two  parallel 

resistances  Ri  and  R?  at  Ti  and  T2  were  measured.  The 

temperature  T2  was  varied  from  293‘’K  to  743®K  while  Ri 

remained  at  room  temperature.  The  linear  dependence  of 

G(f)  on  T2/m  and  the  absolute  value  of  G(f)  were  well 
1 1 


established . 


In  a  second  series  a  capacitance,  C  and  inductance, 

L  were  introduced  in  the  two  branches  of  the  circuit. 
Again  T2  was  varied  and  measurements  made  at  three  dif¬ 
ferent  frequencies  confirming  that  G(f)  depends  only  on 
the  real  part  of  the  branch  impedances  in  accordance  with 
equation  (5) . 

Finally  Ri ,  Rj  and  L  were  kept  at  room  temperature 
and  the  voltage  fluctuations  were  measured  as  a  function 
of  capacitor  temperature.  The  value  of  G(f)  did  not 
change  confirming  that  the  fluctuating  emf ' s  must  be 
associated  only  with  the  resistance. 


3.0  Experimental  Results 

The  primary  objective  of  the  program  is  to  experi¬ 
mentally  verify  equation  (4)  for  solid  and  gaseous  con¬ 
ductors  and  to  demonstrate  that  temperature  can  be  deter¬ 
mined  from  thermal  noise  power  measurements  to  a  working 
accuracy  of  about  20°K  in  the  range  1000°K  to  1600°K. 


3.1  Apparatus 

The  experimental  apparatus  is  shown  schematically  in 
Figure  1 . 

The  high  temperature  environment  is  provided  by  a 
2b“inch  diameter  laboratory  furnace  and  digital  setpoint 
potentiometric  temperature  controller  using  a  chromel- 
alumel  thermocouple.  This  system  can  maintain  a  14  inch 
long  zone  in  the  furnace  to  ±1°K  of  the  indicated  temp¬ 
erature.  The  indicated  temperature  is  monitored  by  direct 
voltage  readout  of  the  thermocouple,  look  up  in  the  manu¬ 
facturers  conversion  table  and  correction  for  room  temp¬ 
erature  variation.  In  this  way  we  estimate  that  the 
temperature  of  the  working  zone  of  the  furnace  can  be 


established  to  within  +2'’K  of  actual. 


An  8  inch  long  1-3/4  inch  diameter  cylindrical  housing 


(Figure  2)  made  of  Inconel  is  used  to  electromagnetically 
isolate  the  hot  resistor  from  external  r.f.  interference. 
The  housing.  Inconel  electrical  leads  and  alumina-silica 
insulating  material  are  all  capable  of  withstanding  the 
maximum  furnace  temperature  of  1600°K. 

In  the  wire  resistor  experiments  a  1  mil  Stablohm 
(70  Ni  30  Cr)  wire  was  attached  from  the  center  lead-in 
to  the  grounded  housing  (Figure  2) .  In  the  gas  resistor 
experiments,  the  air  in  the  annulus  between  the  housing 
and  the  center  lead-in  was  seeded  with  cesium  chloride 
to  provide  a  reistance  somewhere  between  100  and  10,000 
ohms  depending  on  the  temperature  and  vapor  pressure  of 
the  cesium  chloride. 

In  a  1  MHz  bandwidth  the  maximum  available  thermal 
noise  power  ranges  from  4x10  watts  to  just  over  10  ^ 
watts  for  the  300  to  1600°K  range  (Equation  2).  These 
very  small  power  levels  can  be  readily  and  accurately 
measured  using  standard  r.f.  noise  measuring  equipment. 

In  the  present  experiments,  the  noise  from  the  hot  re¬ 
sistor  is  input  into  a  low  noise  amplifier  which  ampli¬ 
fies  the  signal  by  40dB.  The  signal  is  then  passed  onto 
a  spectrum  analyzer  mainframe  which  provides  a  high  qua¬ 
lity  bandpass  filter  centered  at  a  selected  frequency  and 
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then  adds  another  50dB  signal  amplification.  The  output 
from  the  spectrum  analyzer  is  a  21.4  MHz  signal  directly 
pioportional  to  the  filtered  input  power.  This  output 
IS  sent  to  a  thermistor  power  st nsor  which  is  monitored 
by  a  power  meter  with  meter  movement  readout. 

A  v'ector  impedance  meter  with  irn  tor  movomi-nt  readout 
is  used  to  determine  impedance  magnitude  and  phase  for  the 
hot  resistor  experiment  and  the  .imp  i  i  f  ier /power  measuring 
system. 

The  measurement  e  1  ect  i  on  ics  w,  i  .  <  1  b  wed  to  warn  ui 

up  tor  at  least  two  hours  to  ninimw.  d:  itt  in  th.  cali¬ 
bration  settings.  Prior  to  takirn:  ise  powei  nuasuie- 

monts,  the  ca  1  ilirat  ions  ot  th«'  spoct  i  urn  ar.alyzci,  pi  w.  i 
meter  and  vector  impedance  m«  toi  wiio  miielully  chi/CM-d 
and  noted.  We  estimate  that  with  tins  Tuocoduit  iowei 
measurements  are  accurate  to  abeut  1’,  with  most  of  the 
error  attributable  to  reading  thi'  metei  mo\-.-ment. 

One  or  the  largest  sources  ot  r-ot-Mitial  i-rror  in 
these  experiments  is  the  possible  differt*nce  in  tempera¬ 
ture  between  the  location  of  the  t, hei mocoup  1  e  and  the 
interior  of  the  cylindrical  housinu.  To  eliminate  this 
differt'nce  we  allowed  the  cylinder  to  remain  at  tempera¬ 
ture  for  over  30  minutes  until  we  perceived  no  change  in 
measured  power  or  source  impedanci' . 


3.2  Data  Analysis  Techniques 


The  simplest  equivalent  circuit  of  the  measurement 
system  is  shown  in  Figure  3.  In  general  the  impedance 
of  the  temperature  source  is  not  matched  to  that  of 
the  low  noise  amplifier  Zj.  Furthermore  the  thermal 
noise  emitted  by  the  input  of  the  low  noise  amplifier 
must  be  taken  into  account. 

The  response  of  the  amplifier  can  be  represented 
ty  the  set  of  linear  equations 


C  ; 

1  1  1  i 

♦  A  1  ;  1  . 

♦  A  1 

, 

A.  il  , 

♦  A  .  .  1  ; 

+  A 

(7) 

i  V. 

■•xpe 

r  1  Hi’  Tit 

We  Jett* 

rmined 

that 

V;  IS  independent  of 

• 

ind 

t  1.0  t  fi 

e  powei 

sensor 

IS  . 

1  constant  load.  Thus 

A 

-  1) 

and  i 

Thus  equation  (7)  can  be  rewritten 


V  i 

^  ''L 

V 

-  Z.j.i  1 

(8) 

We  consider  that  is  comprised  of  internal  amplifier 

electronic  noise  and  amplifier  input  thermal  noise  vj 

Voc  =  AVl  (9) 
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iFiaure  3)  the  current 
„om  the  general  input  circuit  (Figure 


ii  is 


vs+v_L 

"  2s+Zl 


where  v|  =  AkTgRgAf 


•  ia\  (10) 

Combining  (8) /  '’*  »  ' 


4kAfGTa8^  +  4kAfGT  ^ 

■"  t  ♦»  ^*7  \  *■  ^  S  Lj  * 


l^S+^L' 


whore  Ph  - 


G  = 


11x11 


in  our  experiments  we  measure  P,  Rg  -  ^L' 

.he  valuee  o(  t  anC  hi  are  Known.  The  unKnowne  ^  a^  i-  ^ 
tion  to  the  source  temperature  T,  are  therefore  i. 

.eal  Irts  of  the  source  an.  input  impeaances  an.  are  rn 
general  frequency  dependent. 


V 


To  determine  the  unknowns  we  adopted  the  following 
procedure.  First  we  measured  power  for  several  calibrat¬ 
ed  source  conditions  including  open  circuit,  short  circuit, 
pure  capacitive,  pure  inductive  and  several  calibrated  pure 
resistive  room  temperature  terminations.  These  source  cali¬ 
bration  data  gave  us  reasonably  accurate  values  for  the  4 
amplifier  unknowns. 

We  then  used  equation  (12)  and  the  measured  power  from 
the  hot  resistor  experiments  to  calculate  the  hot  resistor 
temperature,  Tg.  These  calculated  temperatures  are  all 
within  ;.3%  of  the  indicated  tompi-'ratures  measured  by  the 
thermocouple  and  most  of  the  data  is  within  1%.  This  was 
the  highest  accuracy  for  correlating  the  data  that  we  could 
achieve  without  using  a  calibrated  broad  band  noise  source 
with  an  "equivalent  temperature"  of  several  thousand  Kelvins. 

The  most  important  amplifier  parameter  is  the  gain,  G 
since  it  appears  in  two  of  the  terms  of  equation  (12).  The 
most  accurate  value  of  gain  would  best  be  determined  by 
using  an  "equivalent  high  temperature"  calibrated  noise 
source  along  with  our  calibrated  reactive  and  room  tempera¬ 
ture  resistive'  sources. 
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3.3  Final  Results 


The  variation  of  noise  power  with  frequency  bandwidth 
Af  is  shown  in  Figure  4 .  This  data  was  taken  with  wound 
wire  resistors  of  about  160  ohms  for  two  temperatures, 
and  1230°K  and  verifies  the  linear  dependence  of 
noise  power  on  bandwidth. 

The  linear  dependence  of  noise  power  on  resistor  temp¬ 
erature  is  shown  by  the  nickel-chromium  wire  data  in  Figure 
5.  In  these  tests  a  wound  wire  resistor  of  about  75  ohms 
was  used  and  noise  power  was  measured  in  a  3  MHz  bandwidth 
centered  at  20  MHz.  Figure  6  shows  the  calculated  tempera¬ 
ture  based  on  measured  noise  power  versus  the  measured 
thermocouple  temperature.  Figure  7  shows  the  calculated 
noise  power  based  on  measured  thermocouple  temperature 
versus  the  measured  noise  power.  These  correlations  were 
made  using  equation  (12)  .  The  calibration  data  and  ampli¬ 
fier  constants  are  presented  in  Table  1.  The  measured  hot 
wire  data  is  given  in  Table  2. 

The  source  impedance  of  the  hot  wire  experiments  is 
complex  because  of  the  capacitance  and  inductance  of  the 
housing  and  electrical  leads.  However  the  magnitude  and 
phase  of  the  source  impedance  were  constant  to  within  4% 
during  the  experiment. 


Amplif i 


Bandwidth,  Af  [MHz] 


Figure  4.  Noise  Power  from  Hot  Wire  Resistor 
versus  Bandwidth 
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Measured  Temperature  [K] 


Figure  5 


Plot  of  Measured  Amplified 
Noise  Power  versus  Measured 
Hot  Wire  Resistor  Temperature 


Fiqure  6 .  Calculated  Temperature  versus 


Measured  Temperature  for  Hot 


Figure  7.  Plot  of  Calculated  Noise  Power 
versus  Measured  Noise  Power 
Normalized  to  Unit  Bandwidth 
(in  Megahertz) 
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Figure  6.  Plot  of  Calculated  Noise  Power 
versus  Measured  Noise  Power 
Normalized  to  Unit  Bandwidth 
(in  Megahertz) 
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Calibration  Data 


T  =  292°K 


Impedance 

Impedance 

Noise 

Source 

Magnitude 

Phase 

Power 

(ohms) 

(degrees) 

(ywatts ) 

Open  Circuit 

OO 

0 

31.8 

Short  Circuit 

0 

0 

28.4 

Capacitor 

94 

-90 

30.7 

Inductor 

37.9 

90 

31.0 

50  ohm 

50 

.2 

37.4 

75  ohm 

75 

0 

37.2 

600  ohm 

543 

-23.8 

33.4 

Amplifier  Constants 

impedance  Magnitude 

44.2 

ohms 

impedance  Phase 

3.5 

degrees 

Gain 

91.1 

dB 

Frequency 

20 

MHz 

Af 

3 

MHz 

—  Magnitude 

Z  ip 

.01245 

siemens 

4-  Phase 

Zfj* 

183 

degrees 

Pa 

22.02 

ywatts 

Tt 

148.5 

Table  1.  Hot  Wire  Resistor  Calibration 
Data  and  Amplifier  Constants 
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In  the  hot  gas  experiments  however  the  source  imped¬ 
ance  varies  considerably  and  the  data  does  not  correlate 
well  unless  impedance  mismatch  effects  between  source  and 
amplifier  are  accounted  for  properly  as  in  equation  (12). 

The  results  of  the  hot  gas  experiments  are  shown  in 
Figures  6  and  7.  In  these  tests  the  heated  air  seeded 
with  cesium  chloride  formed  the  resistor.  Noise  power 
was  measured  in  a  1  MHz  bandwidth  centered  at  20  MHz.  The 
calibration  data  and  amplifier  constants  are  given  in  Table 
3.  The  measured  hot  gas  data  is  found  in  Table  4. 

The  cesium  chloride  air  resistor  is  a  highly  corro¬ 
sive  agent  on  the  Inconel  surface  of  the  housing  and  elec¬ 
trodes.  In  our  first  experiments  the  surfaces  of  the  elec¬ 
trodes  were  fresh.  In  subsequent  tests  the  walls  were 
coated  with  a  heavy  scale  of  precipitated  cesium  salts. 
However  this  did  not  affect  the  quality  of  the  data.  The 
scale  changes  the  capacitance  of  the  source  but  as  predict¬ 
ed  by  the  theorem  of  Nyquist,  only  the  real  part  of  the 
source  impedance  contributes  to  the  thermal  noise. 
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Calibration  Data 

T  = 

296°K 

Source 

Impedance 

Magnitude 

(ohms) 

Impedance 
Phase 
( Degrees ) 

open  Circuit 

CO 

0 

Short  Circuit 

0 

0 

Capacitor 

92.4 

-90 

I nductor 

37.9 

90 

50  ohm 

50 

.2 

7  5  ohm 

75 

0 

600  ohm 

543 

-23.8 

Noise 

Power 

(vjwatts) 

10.15 

10.4 

9.8 

10.85 

12.6 

12.55 

10.85 


Amplifier  Constants 


impedance  Magnitude 

50.7 

ohms 

Impedance  Phase 

9.6 

degrees 

Gain 

90.71 

dB 

Frequency 

20 

MHz 

af 

1 

MHz 

—  Magnitude 

/jrp 

.00993 

siemens 

A 

y-  Phase 

188 

degrees 

Pa 

7.74 

ywatts 

Tl 

144 

°K 

Table  3.  Hot  Gas  Resistor  Calibration 
Data  and  Amplifier  Constants 
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4.0  Conclusions 


The  primary  objective  of  this  work  is  to  demonstrate 
that  the  temperature  of  hot  gases  such  as  combustion  pro¬ 
ducts  can  be  determined  from  thermal  noise  power  measure¬ 
ments.  This  objective  was  successfully  met  and  we  have 
concluded  that  the  absolute  temperature  can  be  determined 
to  within  1%  when  the  parameters  of  the  amplifying  system 
are  known  to  at  least  this  accuracy. 

This  temperature  sensing  method  can  be  used  to  mea¬ 
sure  quasi-steady  flows  with  response  times  of  about  0.1 
sec  with  currently  available  instrumentation.  Accurate 
submicrosecond  transient  measurements  should  be  possible 
with  an  amplifying  system  based  on  a  Schottky  diode  for 
sensing  noise  power. 
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